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Abstract. This paper presents the modeling of the F9 V star (3 Virginis based on seismological measurements. 
Using the Geneva evolution code including rotation and atomic diffusion, we find that two distinct solutions 
reproduce all existing asteroseismic and non-asteroseismic observational constraints well: a main-sequence model 
with a mass of 1.28 ± 0.03 M0 and an age t = 3.24 ± 0.20 Gyr, or a model in the post-main sequence phase of 
evolution with a lower mass of 1.21 ± 0.02 Mq and an age t = 4.01 ± 0.30 Gyr. The small spacings Svo2 and the 
ratio ro2 between small and large spacings are sensitive to the differences in the structure of the central layers 
between these two solutions and are also sensitive to the structural changes due to the rotational mixing. They 
can therefore be used to unambiguously determine the evolutionary state of /3 Vir and to study the effects of 
rotation on the inner structure of the star. Unfortunately, existing asteroseismic data do not enable such precise 
determination. We also show that the scatter in frequencies introduced by the rotational splittings can account 
for the larger dispersion of the observed large spacings for the non-radial modes than for the radial modes. 
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1. Introduction 

The solar five-minute oscillations have led to a wealth of 
information about the internal structure of the Sun. These 
results stimulated various attempts to detect a similar sig- 
nal on other solar-like stars by photometric or equivalent 
width measurements. In past years, the stabilized spec- 
trographs developed for an extra-solar planet search have 
finally achieved the accuracy needed to make these de- 
tections. While solar-like oscillations have been detected 
for a handful of solar-type stars, individual p-mode fre- 
quencies have only been identified for a few of these stars: 
a Cen A (Bouchy & Carrier f2002: Bedding et al. J004), 
a Cen B (Carrier & Bourban 2003: Kjeldsen et al.|2p05), 
Procyon A (Martic et al. J004a, Eggenberger et al. l2004ai) . 
77 Bootis (Kjeldsen et al. 120031 Carrier et al. I2005a|l . and 
H Arae (Bouchy et al. 12005(1 . Individual frequencies have 
also been identified in the giant star ^ Hydrae (Frandsen 

et ai.nnni- 



Based on these asteroseismic data, numerous theo- 
retical analyses have been performed in order to deter- 
mine precise global stellar parameters (see for instance 
Eggenberger et al. I2004bl for the a Cen system) and to 
try to test the inclusion of overshooting or rotation in the 
models (see Di Mauro et al. i2003i for 77 Boo or Eggenberger 
et al. I2005al for Procyon). 
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Recently, Carrier et al. H2005b|l detected solar-like os- 
cillations on the F9 V star P Virginis (HD 102870) and 
reported identification of several individual frequencies. 
Martic et al. (2004b) also detected solar-like oscillations on 
f3 Virginis, but no individual frequencies were identified. 
The aim of the present paper is to investigate which addi- 
tional constraints are brought by these seismological data. 
We will thus try to determine the model of f3 Virginis that 
best reproduces all the observational constraints using the 
Geneva evolution code, which includes a comprehensive 
treatment of shellular rotation and atomic diffusion. 

The observational constraints available for (3 Virginis 
are summarized in Sect. 2, while the computational 
method is described in Sect. 3. The results are presented in 
Sect. 4 and the effects of rotation are discussed in Sect. 5. 
The conclusion is given in Sect. 6. 



2. Observational constraints 

2.1. Effective temperature and chemical composition 

From the different spectroscopic measurements reported 
in the literature for /3 Vir (Edvardsson et al. 119931 Buzzoni 
et al. nUUT\ Gray et al. l^nUTl Le Borgne et al. HUU^ 
Taylor et al. HUM Morel & Micela HUM AUende Prieto 
et al. 120041 Lambert et al. I2004|l . we adopt an effec- 
tive temperature Tcff = 6130 ± 50K and a metallicity 
[Fe/H] = 0.14 ±0.05. 
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2.2. Luminosity Table 1. Non-asteroseismic observational constraints for 

(i Vir. 

Combining the magnitude V = 3.597 ± 0.004 (Burki et 

al. [2005), the Hipparcos parallax 11 = 91.74 ± 0.77 mas, n [mas] 91 74 zfc 77 

the solar absolute bolometric magnitude Mboi, = 4.746 y [mag] 3.597 ± 0.004 

(Lejeune et al. ll998|l and the bolometric corrections from L/Lq 3.51 ± 0.10 

Flower's l|TM5|l BC = -0.028 ± 0.006 mag determined T^ff [K] 6130 ± 50 

from the effective temperature, we find a luminosity L = [Fe/H] 0.14 ± 0.05 

3.51 ± 0.10 Lq. Fsini [kms-^] ^4.3 



2.3. Rotational velocity 

We use the observed surface velocity of f3 Vir to con- 
strain the rotational velocity of our models. Glebocki & 
Stawikowski H2000|l have determined a velocity Vsini = 
4.3kms~^. Since the value of the angle i is unknown, we 
assume that it is close to 90°. Thus our models of /3 Vir 
have to reproduce a surface velocity of about 4.3 km s""'^. 

All non-asteroseismic observational constraints are 
hsted in Table [2 

2.4. Asteroseismic constraints 

Solar-like oscillations in /3 Vir have been detected by 
Carrier et al. I|20n5b|l with the Coralie echelle spectro- 
graph. Thirty-one oscillation frequencies have been identi- 
fied in the power spectrum between 0.7 and 2.4 mHz with 
amplitudes in the range 23 to 46 cm s~^ . By a least square 
fit of the asymptotic relation with the identified frequen- 
cies, they determined a mean large and small separation 
of Av = 72.1 /L(Hz and 5i/o2 = 6.3 /xHz, respectively. The 
precision of the frequency determination is estimated by 
conducting Monte Carlo simulations in which time series 
computed for various values of the damping time and re- 
excitation rates are analyzed. A velocity time series for de- 
tected oscillation modes is built by using the observational 
time sampling and the observational noise on radial veloc- 
ities. The amplitude spectrum of this series is then calcu- 
lated and the differences between the simulated peaks, and 
the actual frequencies are determined. The whole proce- 
dure is repeated 1000 times for each set of parameters to 
ensure the stability of the results. These simulations lead 
to an estimation of the error on individual frequency be- 
tween 0.8 and 1.4 /xHz, mainly depending on the damping 
time. This result is very close to the frequency resolution 
of the time series (1.14 /itHz). In the following, we thus 
adopt an error on individual frequencies estimated as the 
frequency resolution of the time series. 



3. Stellar models 

3.1. Code description 
3.1.1. General input physics 

The stellar evolution code used for these computations is 
the Geneva code that includes shellular rotation (Meynet 
& MaederEOnnil- We use the OPAL opacities, the NACRE 



nuclear reaction rates (Angulo et al. 11999'), and the stan- 
dard mixing-length formalism for convection. Wc do not 
include an additional convective penetration from the con- 
vective core into the surrounding stable layers, since the 
inclusion of shellular rotation results in an increase of the 
convective core (see Fig. 11 of Carrier et al. I2005a)l . 

3.1.2. Rotation 

We follow the evolution of the rotation profile from the 
zero age main sequence on, assuming initial solid body 
rotation. The braking law of Kawaler (1988) is used in 
order to reproduce the magnetic braking undergone by 
low mass stars when arriving on the main sequence. Two 
parameters enter this braking law: the saturation velocity 
fisat and the braking constant K. Following Bouvier et al. 
1)1997(1 . r^sat is fixed to 141^0 and the braking constant K 
is calibrated on the Sun. 

3.1.3. Meridional circulation 

The velocity of the meridional circulation in the case of 
shellular rotation was initially derived by Zahn H1992)) . 
The effects of the vertical /i-gradient and of the hor- 
izontal turbulence on meridional circulation were then 
taken into account by Maeder & Zahn H1998|l . They find 
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where P is the pressure and Cp the specific heat. The Eq 
and Ef^ terms depend on the f2- and /i-distributions, re- 
spectively, up to the third order derivatives and on various 
thermodynamic quantities (see Maeder & Zahn [1998 for 
more details). 

3.1.4. Shear turbulence 

The diffusion by shear instabilities is expressed by a coef- 
ficient I?shoar: 
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where / is a numerical factor equal to 0.8836, K the ther- 
mal difFusivity, and (V' — V) expresses the difference be- 
tween the internal nonadiabatic gradient and the local gra- 
dient (Maeder & Mevnet .2001f) . 

3.1.5. Horizontal turbulence 

We used the new horizontal turbulence prescription of 
Maeder (2003) that expresses the balance between the 
energy dissipated by the horizontal turbulence and the 
excess of energy present in the differential rotation on an 
equipotential that can be dissipated in a dynamical time: 



with 



Ar {rn{r)V[2V - aU])^ , 
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where V{r) is the vertical component of the meridional 
circulation velocity, V(r) the horizontal component, and 

_ 1 dlnr^O 
" ~ 2 dlnr ' 

3.1.6. Transport of chemicals 

The vertical transport of chemicals through the combined 
action of vertical advection and strong horizontal diffusion 
can be described as a pure diffusive process (Chaboyer & 
Zahn [TM^ . Indeed, the advective transport can be re- 
placed by a diffusive term, with an effective diffusion co- 
efficient 



"soaT 



(5) 



where Z?h is the coefficient of horizontal diffusion (see 
Eq.EJl. 

Atomic diffusion driven by gravitational settling and 
thermal gradients is included using the prescription by 
Paquette et al. (dHEEll- The evolut ion of an element con- 
centration per unit mass Ci then follows: 
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with C/difi the microscopic diffusion velocity and Ci the 
variations of chemical composition due to nuclear reac- 
tions. 

3.2. Computational method 

The characteristics of a stellar model that includes the ef- 
fects of rotation depend on six modeling parameters: the 
mass Al of the star, its age (hereafter i), the mixing- 
length parameter a = I /Hp for convection, the initial 
surface velocity Vi, and two parameters describing the 
initial chemical composition of the star. For these two 



parameters, we choose the initial hydrogen abundance 
Xi and the initial ratio between the mass fraction of 
heavy elements and hydrogen {Z/X)i. This ratio is di- 
rectly related to the abundance ratio [Fe/H] assuming that 
log{Z/X) = [Fc/H] + \og{Z/X)Q; we use the solar value 
{Z/X)rs, = 0.0230 given by Grevesse & Sauval ifT^ . 

Contrary to our previous seismic calibrations of bi- 
nary stars like a Centauri and Procyon, we decided to 
limit the parameter space of the present analysis by fixing 
the mixing-length parameter to its solar calibrated value 
(a0 = 1.75) and the hydrogen abundance to X — 0.7. 
This simply results from the limited number of non- 
asteroseismic observational constraints available for a star 
like /3 Virginis, which does not belong to a binary system. 
Although this arbitrarily limits our analysis to a small sub- 
set of possible solutions, we think that it is the appropri- 
ate way to proceed for a first modeling of an isolated star 
in order to investigate which constraints are brought by 
the observed p-mode frequencies. A more thorough explo- 
ration of the parameter space, leading to the independent 
determination of the mixing-length parameter and the hy- 
drogen abundance, will be possible when additional and 
more accurate seismic and non-asteroseismic observables 
(e.g. interferometric radius) are available for /3 Virginis. 

We then construct a grid of models with position in the 
HR diagram in agreement with the observational values of 
the luminosity and effective temperature (see Fig.QJ. Note 
that the initial ratio between the mass fraction of heavy el- 
ements and hydrogen {Z/X)i is directly constrained by the 
observed surface metallicity [Fe/H], while the initial veloc- 
ity Vi is directly constrained by the observed rotational ve- 
locity. For each stellar model of this grid, low-£ p-mode fre- 
quencies are calculated using the Aarhus adiabatic pulsa- 
tions package written by J. Christensen-Dalsgaard H1997|l . 
Following the observations, only modes £ < 2 with fre- 
quencies between 0.7 and 2.4 mHz are computed. To de- 
termine the set of modeling parameters {M, t, Vi, {Z/X)i) 
that leads to the best agreement with the observational 
constraints, two functionals are defined: Xcias ^^'-^ xLtcro- 
The Xcias functional only uses the "classical" (i.e. non- 
asteroseismic) observables, and is defined as follows 



Xcla 



-j^ 3 / ^tlieo (^obs ^ ^ 



i=l 



obs 
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where the vector C contains the following observables for 
one star: 

C = (i/L0,Teff,[Fe/H]) . 

The vector 0**^°° contains the theoretical values of these 
observables for the model to be tested, while the values of 
Qobs ^]-^Qgg listed in Table The vector trC contains 
the errors on these observations. Note that the observed 
rotational velocity is not included in this minimization, be- 
cause of its large uncertainty resulting from the unknown 
inclination angle i. The xLtcro functional compares indi- 
vidual theoretical frequencies to the observed ones. The 
exact values of the computed p-mode frequencies depend 
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on the details of the star's atmosphere, where the pulsa- 
tion is non-adiabatic and where turbulent and radiative 
losses are significant. The discrepancy between observed 
and theoretical solar p-mode frequencies calculated from 
standard models is indeed caused primarily by the approx- 
imations made in modeling the surface layers (Ulrich & 
Rhodes il983l Christensen-Dalsgaard et al. ll988l Guenther 
n^MI Li et al. EnnUl- The effects of poor modehng of the 
external layers on computed p-mode frequencies can be 
taken into account by using the discrepancies between ob- 
served and theoretical solar frequencies (see Guenther & 
Brown|^04). However, we do not know whether this esti- 
mation is reliable for stars with different surface gravity, 
chemical composition, and age. Thus, we prefer to sim- 
ply apply a constant shift to correct theoretical frequen- 
cies for surface effects. We are well aware that this simple 
correction may seem doubtful, given the strong frequency 
dependence of the solar near-surface effects, but we think 
that this is the least biased procedure possible. We there- 
fore define the mean value of the differences between the 
theoretical and observed frequencies: 



if r 



thco 



,obs 



i=l 



and then define the xLtcro functional as 
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where a = 1.14 /xHz is the error on the observed frequen- 
cies, estimated as the frequency resolution of the time 
series (see Sect. 12. 4|) . and N is the number of observed 
frequencies. Note that the adopted value of the error on 
the observed frequencies will of course change the values 
of the computed xLtcro larger value of u resulting in 
smaller values of xLtcro ^iid vice versa) but not the result 
of the minimization of the xLtero functional. Carrier et al. 
Ij2005b|l report the identification of thirty-one individual 
frequencies. However, four modes are identified at high 
frequencies and are thus far from the other modes, which 
makes their identification less secure. Moreover, these fre- 
quencies are found to be larger than the typical value of 
the acoustical cut-off frequency of our models of j3 Vir 
(z/ac — 2.1 mHz). As a result, these four frequencies are 
not considered in our minimization, and N is hence equal 
to 27 instead of 31. 

Determination of the best set of parameters is then 
based on the simultaneous minimization of the two func- 
tional defined in Eqs. [7| and|Sl Finally, the global agree- 
ment between the stellar models and all classical and as- 
teroseismic observational constraints available for (3 Vir is 
tested by defining the Xtot functional as the sum of x, 

Xastcro- 



clas 



2 — 2 
Xtot Xclas 



Aastcro 



(9) 
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Fig. 1. Evolutionary tracks in the HR diagram for the two 
models of j3 Vir. The dot and the square indicate the loca- 
tion of the Ml and M2 models, respectively. The error box 
in continuous line indicates the observational constraints 
of L and Toff. The evolutionary tracks of a 1.19 Af© and 
1.31 Mq model are also plotted (dotted lines); this shows 
that models with masses that are lower than 1.19 Af© or 
larger than 1.31 Af© are not able to reproduce the non- 
asteroseismic constraints. 



4. Results 

Using the observational constraints listed in Sect. |21 with 
the observed frequencies of Carrier et al. H2005b|l . we per- 
form the x^ minimization described above. 

First, only the non-asteroseismic constraints are con- 
sidered. Thus the Xcias minimization, which tests the con- 
sistency between observed and theoretical effective tem- 
perature, luminosity and metallicity, is applied. In this 
way, we find that two kind of solutions are able to repro- 
duce the non-asteroseismic observational constraints of (3 
Vir: one model with a mass of about 1.21 Mq and an age of 
about 4.1 Gyr, and a more massive model of 1.28 Af© with 
a smaller age of about 3.2 Gyr. This is shown in Fig. |21 
which plots the values of x?ias ^o^' models matching the 
observed surface metallicity of (3 Vir as a function of age 
and mass. This figure clearly shows one minimum near 
1.21 Af© (corresponding to the minimum at 4.1 Gyr) and 
another minimum near 1.28 A/© (3.2 Gyr). To understand 
why two distinct solutions arc found to fall within the 
observed constraints on effective temperature, luminosity 
and surface metallicity, we plot the evolutionary tracks in 
the HR diagram corresponding to the 1.21 and 1.28 Ai"© 
models (see Fig.^. Figure^shows that the two solutions 
correspond in fact to two different evolutionary stages: the 
less massive model of 1.21 Af© has exhausted its hydro- 
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gen core and is therefore in the post-main sequence phase 
of evolution, while the 1.28 Mq is still on the main se- 
quence but is very near hydrogen exhaustion. Since these 
two models have the same surface metallicity, the same 
surface velocity {V ^ 4knis^^), and the same position 
in the HR diagram, non-asteroseismic observational con- 
straints do not enable us to determine the evolutionary 
stage of /3 Vir. Asteroseismic observations are needed to 
differentiate models with different internal structure lo- 
cated in the same region of the HR diagram. Thus, we 
now examine which additional constraints are brought by 
the asteroseismic measurements by performing the whole 
minimization described in Sect.|2| 




age [Gyr] mass [Mg] 



Fig. 2. Top: x^ias ^^i*^ models as a function of age and 
mass. Bottom: Zoom of the above figures. Models with 
Xastero > Still denoted by crosses, while models with 

Xastero — 10 indicated by dots. 

For this purpose, the Xastoro values, which directly 
compare individual theoretical frequencies to the observed 
ones, are computed. The bottom of Fig. |21 shows which 
models are in reasonable agreement with the observed p- 
mode frequencies (with xLtoro ^ 10) ^^'^ which models do 
not reproduce the asteroseismic constraints (xLtcro > 10)- 
We see that models still on the main-sequence, which min- 
imize Xcias ^'^^ ^ mass and a corresponding age of about 
1.28 and 3.2 Gyr respectively, are found to also agree 
well with the seismological measurements, and hence to 
minimize x?ias ^^'^ Xastero the same time. However, 
models in the post-main sequence phase of evolution, 
which minimizes Xcias ^'^^ ^ mass and a corresponding age 
of about 1.21 M0 and 4.1 Gyr, have values of xLtcro ^^^^ 
are higher than ten and therefore badly reproduce the as- 
teroseismic observables. Thus, post-main sequence mod- 



els are not able to minimize the two functionals Xcias ^^^'^ 
Xastero ^.t the Same time as expected for a fully consistent 
model of (3 Vir. However, the bottom of Fig.Elshows that 
there are post-main sequence models that are compatible 
with the observed frequencies and that also have small 
values of Xcias' even if they do not minimize Xcias- This 
means that a post-main sequence model of /3 Vir exists 
with position in the HR diagram in agreement with the 
observed luminosity and effective temperature, which also 
reproduces the observed p-mode frequencies well. This can 
be seen clearly by plotting the lowest values of the Xtot 
functional, which tests the global agreement between the 
models and all observational constraints, as a function of 
age and mass (see Fig. EJ. Figure 01 shows that there is 
a clear solution near 1.28 M© that minimizes all x^ func- 
tionals at the same time. Figure 01 also shows that the 
post-main sequence models near 1.21 Mq have low values 
of Xtot ^n.d are therefore in good agreement with all ob- 
servational constraints. At this stage of the analysis, we 
conclude that the additional asteroseismic measurements 
favor the solution of a model still on the main sequence, 
but do not rule out the other solution of a less massive 
model in the post-main sequence phase of evolution. 




age [Gyr] mass [M^] 



Fig. 3. Lowest values of xL = Xcias + xLtero as a function 
of age and mass. Dots indicate models in accordance with 
the observed luminosity and effective temperature, while 
open circles correspond to models that do not reproduce 
all non-asteroseismic observables. 

Thus, two different solutions are determined: a main- 
sequence model (denoted model Ml in the following) and 
a post-main sequence model (the M2 model). For each 
model, the location in the HR diagram is shown in Fig.^ 
and the respective characteristics given in Table [3 The 
confidence limits of each modeling parameter are esti- 
mated as the maximum/minimum values that fit the ob- 
servational constraints when the other calibration param- 
eters are fixed to their mean values. Both models have a 
similar value of Xastero i while the M2 model has a higher 
value of Xcias' which is mainly due to its high effective 
temperature (see Fig^. 

We recall that these results have been obtained by fix- 
ing the mixing-length parameter to its solar calibrated 
value (a© — 1.75) and the hydrogen abundance to X — 
0.7, since we think that current observational data do not 
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Table 2. Models for (3 Vir including rotation and atomic 
diffusion. The upper part of the table gives the non- 
asteroseismic observational constraints used for the cal- 
ibration. The middle part of the table presents the mod- 
eling parameters with their confidence limits, while the 
bottom part presents the global parameters of the star. 





Model Ml 


Model M2 


T 1 T 

L/Lq 


3.51 ±0.10 


Teff [K] 


6130 ± 50 


[Fe/H] 


0.14 ±0.05 


V [kms~-^] 




4.3 


M [Mq] 


1.28 ±0.03 


1.21 ±0.02 


t [GyA 


3.24 ±0.20 


4.01 ±0.30 


Vi [kins ^] 


~ 18 


~ 18 


(Z/X)i 


0.0340 ± 0.0040 


0.0340 ± 0.0040 


Xclas 


0.002 


0.319 


2 

Xastcro 


1.961 


2.003 


2 

Xtot 


1.963 


2.322 


L/Lq 


3.51 


3.49 


Teff [K] 


6126 


6177 


R/ i?Q 


1.666 


1.634 


[Fe/H] 


0.14 


0.14 


(Z/X) 


0.0317 


0.0317 


V [kms"'] 


4.2 


3.9 


Av [/iHz] 


72.1 


72.1 


6vo2 [mHz] 


5.4 


5.5 



allow an independent determination of these quantities. 
We nevertheless investigated the effect of a change to these 
values on the determined global parameters. Increasing 
the mixing-length parameter a by 10 % for a given value 
of the hydrogen abundance results in an increase of about 
0.1 in the x^ias order to reach the same value of xLtoroi 
this also decreases the mass by about 2 % and increases 
the age by about 10 %. In the same way, increasing the hy- 
drogen abundance by 10 % for a fixed value of the mixing- 
length parameter leads to an increase of about 0.7 in the 
Xclas order to reach the same value of xLtcro; this in- 
creases the mass by about 12 % and decreases the age by 
about 6%. For a given mass, an increase in the hydro- 
gen abundance can also be compensated for an increase 
in the mixing-length parameter and vice versa in order to 
reach the same location in the HR diagram (see Fig. 3 of 
Eggenberger et al. |2005a). 

To study the asteroseismic features of the models, we 
first plot the variation of the large spacing {Aunj = 
Vn,t — i^n-i,t) with frequency; the comparison between ob- 
served and theoretical large spacings for the Ml and M2 
model is given in Fig. 0] Table |21 and Fig. 0] show that 
the mean large spacing of both models are in good agree- 
ment with the observed value of 72.1 /iHz. Note that the 
radius of the M2 model is less than the one of the Ml 
model in order to compensate for the decrease in mass 
and to reproduce the observed value of the mean large 
spacing. The variations in the theoretical large spacings 
with frequency are also very similar for both models. The 
only difference is that the M2 model exhibits slightly lower 



values of the £ = 1 large spacings at low frequency than 
does the Ml model. As a result, the Ml model is able to 
reproduce the observed t — 1 large spacing at 800 /iHz, 
whereas the theoretical value of the M2 model is lower 
than the observed one. This is due to the fact that the 
M2 model does not reproduce the observed t — 1 mode at 
732 /iHz well (see Fig.[71l. Finally, we note that the disper- 
sion of the observed large spacings around the theoretical 
curves is greater than expected taking an uncertainty of 
1.14yLtHz, the frequency resolution of the time series, on 
the frequency determination. 




Fig. 4. Large spacings versus frequency for the Ml and 
M2 models. Dots indicate the observed values of the large 
spacing for radial modes, while triangles and squares cor- 
respond to the observed large spacings determined from 
i —\ and i — 2 modes, respectively. 

The small spacing 5vo2 = i'n.i=a — i^n-i.e=2 is prin- 
cipally sensitive to the conditions in the central regions 
of the star. However, it also retains some sensitivity to 
the near-surface structure. Roxburgh & Vorontsov H2()03|l 
introduces the use of a new asteroseismic diagnostic: the 
ratio of small to large spacings ro2 = 5i^02/At'„/=i. This 
ratio constitutes a better diagnostic of the central parts 
of a star than the small spacing, since it is essentially 
independent of the structure of the outer layers, and is 
determined solely by the interior structure (Roxburgh & 
Vorontsov EOOSI Oti Floranes et al. lMT^ . 

The variation in the small spacing Siyo2 and the ratio 
ro2 with frequency for the Ml and M2 models is given in 
Fig. 13 The large spacing corresponding to the observed 
small spacing near 1700 /iHz has not been observed; thus, 
only two observational values of the ro2 ratio are available. 
We see that the theoretical mean small spacing of both 
models (see Table EJ is smaller than the observed value of 
6.3±1.4 /xHz. Note that the theoretical mean small spacing 
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ration ro2 = 5vq2/ ^Vn^i^i as a function of frequency. The 
continuous line corresponds to the Ml model, while the 
dotted line corresponds to the M2 model. Dots indicate 
the observed values. 



has been determined by only considering the three theoret- 
ical values corresponding to the observed small spacings, 
in order to correctly compare stellar models to observa- 
tions. Figure El shows that the difference between the the- 
oretical and observed mean small spacings is mainly due to 
the large observational value of 5vq2 near 1270 /iHz, which 
models fail to reproduce. This is also true for the ro2 ratio. 
The variation in Svm and ro2 with frequency is very sim- 
ilar for both models at low frequency, but is found to be 
quite different at higher frequency. Indeed, Fig. [S] shows 
that the Ml and M2 models are characterized by approxi- 
mately the same values of 51/02 and ro2 at frequencies lower 
than about 1600 /itHz, while at higher frequencies, the Ml 
model exhibits a significantly stronger decrease in (5i'o2 
and ro2 than does the M2 model. This is directly related 
to the evolutionary stage of the two models. Indeed, the 
post-main-sequence M2 model is more centrally condensed 
and therefore exhibits higher central values of pressure and 
density than the younger Ml model, which is approach- 
ing hydrogen exhaustion but is still on the main-sequence. 
Consequently, the sound speed and hence the Lamb fre- 
quency of the two models differ in the central parts of the 
star. The interior reflection of the 1 = 2 modes are then 
different for the two models, which results in a stronger de- 
crease in 5vo2 and ro2 for the Ml model than for the M2 
model. Unfortunately, only one value for 5vq2 and none 
for ro2 have been observed for a frequency higher than 
1600 /iHz. The Ml model agrees slightly better with this 
value of the small spacing near 1700 /iHz than does the 
M2 model, but the difference is clearly not significant. 



Finally, theoretical and observed p-mode frequencies 
are compared by plotting the echelle diagram of the Ml 
and M2 models (Figs. El and d). In these figures the sys- 
tematic difference {D^) between theoretical and observed 
frequencies has been taken into account {{Di,) = 23.8 /iHz 
and 21.1 /xHz for the Ml and M2 model, respectively). 
This shift of about 20 /iHz between theoretical and ob- 
served frequencies is substantially larger than the shift 
required in the solar case. We recall that our models have 
been computed with the standard mixing-length formal- 
ism for convection and that we assume the oscillations to 
be adiabatic. It is beyond the scope of the present paper 
to investigate the influence of these simple assumptions on 
the solutions, but we note that they may explain a shift of 
about 20 /iHz between computed and observed frequencies 
for a star that is more massive and more evolved than the 
Sun. 
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Fig. 6. Echelle diagram for the Ml model. Open symbols 
refer to the theoretical frequencies, and filled symbols to 
the observed ones. Circles are used for ^ = modes, tri- 
angles for t = 1 modes, and squares ior i = 2 modes. The 
two £ = 2 modes near 1480 /iHz and separated by 1.4 /iHz 
could be due to rotational splitting. 



5. Effects of rotation 

5.1. Small spacings 5vq2 anc/ ratio ro2 

As already mentioned, the theoretical mean small spacings 
of the Ml and M2 models are smaller than the observed 
ones. The theoretical values are, however, perfectly com- 
patible with the observations given the large uncertainty 
on the observed mean large spacing. Moreover, only three 
values of 5vq2 (and two of ro2) have been observed. It is 
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nevertheless interesting to investigate which physical pro- 
cess can increase the theoretical value of the mean small 
spacing in order to improve the agreement between the 
models and the seismological data. 

Because the small spacing is principally sensitive to 
the structure of the core, its mean value decreases dur- 
ing the evolution of the star on the main sequence. Thus, 
younger models than the Ml model, which is very near 
hydrogen exhaustion, exhibit higher values for the mean 
small spacing. Although these younger models are in bet- 
ter agreement with the observed values of Si'q2 , they fail to 
reproduce the other seismological data. Indeed, the mean 
large spacing also decreases during the evolution of the 
star on the main sequence due to the increase of its ra- 
dius. Consequently, the mean large spacing of the younger 
models is larger than the observed one; these models are 
therefore unable to reproduce the observed individual fre- 
quencies. 

It is thus necessary to find a process that is able to in- 
crease the mean small spacing without changing the mean 
large spacing. By changing the conditions in the central 
regions of the star, the rotational mixing is able to in- 
crease the values of 6i'q2 and ro2 for models with the same 
mean large spacing. To investigate the effects of rotational 
mixing on Si'02 and r'02, we compute models with exactly 
the same initial parameters as the Ml model, but differ- 
ent initial velocities. The Ml model is computed with an 
initial velocity of 18kms^^ in order to reproduce the ob- 
served surface velocity of about 4 km s^ ^ . We calculate two 
other models with higher initial rotation velocities of 50 
and lOOkms"^. Both models are evolved until their mean 



large spacing equals the observed value of 72.1 /iHz. The 
surface velocity V decreases during the evolution on the 
main sequence as a result of the magnetic braking that 
undergoes the low mass stars. When the observed value of 
the large spacing is reached, these velocities are equal to 
6.8 and 8.2 km s"^ for the model with an initial velocity 
of 50 and 100 km s^^, respectively. 

In this way, we obtain three models computed with the 
same input physics and initial parameters except for the 
initial velocity, and they are in good global agreement with 
the asteroseismic measurements. The models with a higher 
initial velocity are of course in poorer agreement with the 
classical constraints than the Ml model, since an increase 
in the rotational velocity results in an increase in the effec- 
tive temperature and luminosity. All models are, however, 
compatible with the adopted non-asteroseismic observa- 
tional constraints. The variation in the small spacing Si'02 
and the ratio ro2 with frequency for the three models is 
given in Fig.|Hl This figure shows that an increase in the ro- 
tational velocity increases di>02 and ro2. Indeed, the mean 
small spacing of the Ml model is 5.4 /iHz, while it is equal 
to 5.7 and 6.0 /iHz for the model with a surface velocity of 
6.8 and 8.2kms~^, respectively. This illustrates the fact 
that the mixing in the central layers, hence the values of 
(5i^02 and ro2, increase with the rotational velocity. 

5.2. Surface abundances 

In addition to changing the structure of the stellar core, 
rotation also influences the chemical abundances at the 
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surface of the star. Indeed, rotation-induced mixing coun- 
teracts the effects of atomic diffusion in the external lay- 
ers of the star. Since the efficiency of this mixing increases 
with the rotational velocity, the decrease of the helium 
and heavy elements abundances due to atomic diffusion is 
found to be smaller for larger rotational velocities. This 
is illustrated in Fig. |51 which shows the helium profile in 
the external layers of the star for the three rotating mod- 
els with the same mean large spacing of 72.1 /iHz. We see 
that the helium surface abundance is larger when the sur- 
face velocity is higher. Although the effects of rotation on 
the surface abundances are relatively small for stars with 
masses lower than 1.3 like fi Vir, we note that they 
are particularly important for more massive stars with 
shallower convective envelope like Procyon or 77 Bootis. 
Indeed, in stars that are more massive than about 1.4 Mq, 
the rotation induced mixing prevents the helium and the 
heavy elements from being drained out of the convective 
envelope (see Fig. 10 of Carrier et al. l2()()5al for rj Bootis). 
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Fig. 9. Helium abundance profile in the external layers of 
the star for models with the same value of the mean large 
spacing (72.1 ^Hz) but different surface velocities V. 
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Fig. 11. Rotational splittings {5vrot{n,(.) = 
Sujn,e,m=i/'2T^) for the Ml model. The dotted and 
the dashed lines correspond to the splittings for modes 
with £ ^ 1 and £ = 2, respectively. 



5.3. Dispersion of the observed large spacings 

As mentioned above, the observed large spacings show 
a significant dispersion around the theoretical curves. 
Carrier et al. ((2005bfl suggest that this scatter can be 
partly explained by rotational splitting, since the disper- 
sion for £ = 2 and £ — 1 modes are found to be greater 
than for radial modes. To check this, we compute the the- 
oretical rotational splittings for the £ — 2 and £ = 1 modes 
of the Ml model. These frequency splittings are calculated 
using 

SuJnim^mPne Kni{r)Q.{r)dr , (10) 
Jo 

where the rotational kernel Kni and the coefficients /?„£ 
are defined by Eqs. (8.42) and (8.43) from Christensen- 
Dalsgaard H2003II . The rotation profile of the Ml model 
and an example of the rotational kernel computed for the 
£ = 2 mode with n = 16 are shown in Fig. ^| Using Eq. 
Hl()(l . the rotational splittings (Ji'rot expressed in terms of 
cyclic frequencies {5v-cot{n,£) = (5w„_£_m=i/27r) are calcu- 



The theoretical rotational splittings lie between 0.8 
and 0.9 /iHz for ^ = 1 modes, and between 0.75 and 
0.85 /iHz for £ — 2 modes. We notice that these values 
are higher than the value of about 0.6 /iHz deduced by as- 
suming a uniform rotation velocity equal to the angular 
velocity at the surface. This reflects the increase of the 
angular velocity SI when the distance to the center of the 
star decreases (see Fig. IIOII . Due to the rotational split- 
tings, the frequencies of the t = 1 modes can be shifted 
by ±5v^ot{£ — 1) (~ 0.85 /iHz), and the frequencies of the 
i — 2 modes can be shifted by ±25v-[ot{£ — 2) 1.6 /iHz). 
This introduces a scatter of the frequencies which can ex- 
plain the high dispersion of the observed large spacings. 
Figure E| shows the theoretical dispersion of the large 
spacings expected from rotational splittings; the values for 
the £ = \ modes are delimited by the modes with m = ±1, 
while the values for the t — 2 modes are delimited by the 
modes with m — ±2. 

We see that the scatter of frequencies resulting from 
the rotational splittings can account for the high disper- 
sion of the observed large spacings and especially for the 
fact that this dispersion is greater for £ = 2 and £ ~ 1 
modes than for radial modes. Indeed, the observed large 
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Fig. 12. Dispersion of the large and small spacings ex- 
pected from rotational splittings for the Ml model. The 
observed large spacings for the £ = 1 modes are expected 
to lie between the maximum/minimum values delimited 
by the I — \^ m — ±1 modes. In the same way, the ob- 
served large spacings for the 1=2 modes and the observed 
small spacings &vq2 are expected to lie between the max- 
imum/minimum values delimited by the i = 2, m = ±2 
modes. 



spacings for the 1 = 1 and 1 = 2 modes are found to 
lie within the predicted interval in Ai^. There is only 
one exception: the i = 1 value near 1880 //Hz, which is 
too small to be correctly reproduced by the models. We 
also notice that the two values of the i = large spac- 
ing near 1400 /iHz are not in good agreement with the 
theoretical predictions. This cannot, of course, been ex- 
plained by rotational splittings, but can either be due to 
the scatter resulting from the finite lifetimes of the modes 
(which, contrary to the rotational splittings, affects all the 
modes including the radial ones) or/and to the misidenti- 
fication of the observed frequency near 1340 /iHz. Indeed, 
the echelle diagrams show that the models fail to repro- 
duce the (. = frequency near 1340 /iHz (see Figs. Eland 
d). As a result, the two observed large spacings for I = Q 
modes near 1400 /iHz do not agree well with the theoretical 
ones. It is interesting to note that the extracted frequency 
corresponding to the ^ = frequency at 1340.7/iHz is 
1329.1 /iHz (see Table 2 in Carrier et al. l2005b|l . Thus, this 
frequency has been shifted by -)-11.57/iHz by Carrier et 
al. and has been identified as an £ = mode. Comparison 
with our models shows that the extracted frequency of 
1329.1 /iHz in fact better agrees with the theoretical pre- 
dictions than does the shifted frequency of 1340.7/iHz. 
Thus, we think that this frequency, which was originally 
identified as a radial mode at 1340.7/iHz by Carrier et al., 



is in fact an £ = 2 mode at 1329.1 /iHz. The observed small 
spacings are also perfectly compatible with the theoretical 
predictions, when the uncertainties due to the rotational 
splittings of the t = 2 modes are taken into account (see 
bottom of Fig. 112). 

Figure^! shows that models that include shellular ro- 
tation predict an increase in the angular velocity when the 
distance to the center decreases in contrast to the near 
uniformity of the solar rotation profile. Indeed, helioseis- 
mological results indicate that the angular velocity VL{r) 
is constant as a function of the radius r between about 
20 % and 70 % of the total solar radius (Brown et al. [T!M?l 
Kosovichev et al. 119971 Couvidat et al. 2003). Since shel- 
lular rotation alone produces an insufficient internal cou- 
pling to ensure solid body rotation (Pinsonneault et al. 
1989: Chaboyer et al. lT^ Eggenberger et al. l2005b|l . this 
suggests that another effect intervenes. Two mechanisms 
have been proposed as reproducing the solar rotation pro- 
file: magnetic fields (Mestel & Weiss [TpFf Charbonneau 
& MacGregor .1993., Eggenberger et al. Il2005bj and inter- 
nal gravity waves (Zahn et al. 119971 Talon & Charbonnel 
P005). However, the question remains open as to how ro- 
tation should be modeled in a slightly more massive solar- 
type star like (3 Virginis. It is therefore worth investi- 
gating how a flat rotation profile will affect our conclu- 
sions about the effects of rotation on the small separa- 
tions and the rotational splittings of [3 Virginis. The near 
solid body rotation of models including magnetics fields 
results in a low value for the diffusion coefficient asso- 
ciated to the shear turbulent mixing, and suggests that, 
for slow rotating solar- like stars, rotational mixing is less 
efhcient in magnetic models than in models with only ro- 
tation (Eggenberger et al. I2005b| . In the same way, the 
inclusion of internal gravity waves is found to reduce ro- 
tational mixing together with differential rotation (Talon 
& Charbonnel I2005|l . The mixing in the central layers for 
a model of (3 Virginis exhibiting a solar-like rotation pro- 
file will thus be reduced resulting in a less pronounced 
increase in 5va2 and ro2 with the rotational velocity. As 
mentioned above, the values of the rotational splittings 
will also be influenced: the rotational splittings computed 
for a model of (3 Virginis with a uniform rotation veloc- 
ity are smaller than the splittings corresponding to the 
rotation profile of our Ml model (mean value of 0.57 /iHz 
instead of 0.83 /iHz). This results in a smaller scatter in the 
frequencies of non-radial modes for models with a solar- 
like rotation profile than for our Ml model. It is, how- 
ever, clear that the accuracy of the present asteroseismic 
observations does not enable us to distinguish between 
these models with different rotation profiles. We can sim- 
ply hope that future asteroseismic observations will lead to 
accurate determination of rotational splittings and small 
separations for stars with various masses and ages, in or- 
der to provide us with new insight into the rotation of 
solar-type stars and with constraints on how to model it 
best. 
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6. Conclusion 

The aim of this work was to determine the best model for 
the F9 V star (3 Virginis using the Geneva evolution code. 
By combining the existing non-asteroseismic observables 
with the new seismological measurements by Carrier et 
al. (2003B, we find that two different solutions reproduce 
all these constraints well: a main-sequence model with 
a mass of 1.28 ± 0.03 Mq, an age t = 3.24 ± 0.20 Gyr, 
and an initial metaUicity {Z/X)i = 0.0340 ± 0.0040, or 
else a model in the post-main sequence phase of evolution 
with a lower mass of 1.21 ± 0.02 and a larger age t = 
4.01 ± 0.30 Gyr. Although these two models are compati- 
ble with all observables, we note that the main-sequence 
model agrees slightly better with the non-asteroseismic 
observables than the post-main-sequence one. 

We show that the small spacings 5vo2 and the ratio ro2 
between small and large spacings are sensitive to the dif- 
ferences in the structure of the central layers between the 
main-sequence and the post-main-sequence model; there- 
fore, they can be used to unambiguously determine the 
evolutionary stage of /3 Vir. Unfortunately, existing aster- 
oseismic data do not enable such a determination, since 
none of the observed 6vq2 and ro2 are available for fre- 
quencies where significant differences between the values 
of (5i/o2 and ro2 for main-sequence and post-main-sequence 
models are predicted. 

The effects of rotation were also studied. We first find 
that the changes in the inner structure of the star due to 
the rotational mixing may be revealed by asteroseismic 
observations, if precise measurements of 5vo2 and ro2 can 
be made. We also show that the scatter of frequencies 
introduced by the rotational splittings can account for the 
larger dispersion of the observed large spacings for the 
non-radial modes than for the radial modes 

We conclude that the combination of non-asteroseismic 
measurements with existing frequency observations puts 
important constraints on the global parameters of (3 Vir, 
but that additional asteroseismic measurements are 
needed (especially measurements of the small spacings) 
to unambiguously determine its evolutionary stage and to 
investigate the internal structure of this star in detail. 
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